Quantum magneto-optics of graphite with trigonal warping 
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The optical conductivity of graphite in quantizing magnetic fields is studied. Both the dynamical 
conductivities, longitudinal as well as Hall's, are analytically evaluated. The conductivity peaks are 
explained in terms of electron transitions. We have shown that the trigonal warping in graphite 
can be considered within the perturbation theory at the strong magnetic field larger than 1 T 
approximately. The main optical transitions obey the selection rule with An = 1 for the Landau 
number n, however the An = 2 transitions due to the trigonal warping with the small probability 
are also essential. The Kerr rotation and reflectivity in graphite in the quantizing magnetic fields 
are calculated. Parameters of the Slonczewski-Weiss-McClure model are used in the fit taking into 
account the previous dHvA measurements and correcting some of them for the case of the strong 
magnetic fields. 

PACS numbers: 76.40.+b, 78.20.-e, 81.05.uf 
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I. INTRODUCTION 

Properties of graphite have attracted much atten- 
tion for more than 50 years. Many of that proper- 
ties were successfully explained within the Slonczewski- 
Weiss-McClure (SWMC) theory^. The most accurate 
method to study the band structure of graphite is a 
study of the Landau levels (LLs) through experiments 
such as magneto-optics 2 - - — and magneto-transport^— . 
However, the interpretation of the experimental results 
involves a significant degree of uncertainty since, as it 
is not clear how the resonances should be identified and 
which electron transitions they correspond to. 

The SWMC theory requires the use of many tight- 
binding parameters and provides the simple description 
of observed phenomena either in the semiclassical limit 
of week magnetic fields or for high frequencies when the 
largest tight-binding inter-layer parameter 71 plays the 
leading rolei^. It is more difficult to take into account the 
parameter 73 known as " trigonal warping" . Usually, it 
is either neglectec&i£r— or considered numerically 2 ^—. 
The Bohr-Sommerfeld quantization condition was ap- 
plied in ReL2£ to find in low magnetic fields the level 
structure including the trigonal warping. In any case, 
only the problem of levels was considered so far, and 
no calculations of conductivities were done in order to 
evaluate the optic properties of graphite. The problem 
appearing for three-dimensional systems in the magnetic 
field connects partly with integrating over the momen- 
tum projection k z along the magnetic field. 

The SWMC model can be simplified assuming that 
only the integration limits such as the K and H points in 
the Brillouin zone produce the main contributions 8 - : 1 ' 18 . 
Such an approximation is similar to the theory of 
magneto-optical effects in topological insulators 2 ^ and 
graphene 2 ^. However, the band extrema or the integra- 
tion limits at the Fermi level can also contribute into the 
absorption. Therefore, the analytical expression for the 
dynamic conductivity in the presence of magnetic fields 



is needed for the interpretation of the magneto-optics ex- 
periments. 

In this paper, motivated by the experimental study 
of the Faraday rotation in single- and multilayer 
graphenei^, we propose a theory of magneto-optics phe- 
nomena in graphite in strong magnetic fields including 
the interlayer hopping parameters 73 and 74 in the Hamil- 
tonian. The trigonal warping 73 is considered as a per- 
turbation with the help of the Green's function method. 
Not only the energy-level structure corrected due to the 
trigonal warping is found, but the expressions for the lon- 
gitudinal and Hall dynamical conductivities are derived. 
Our main theoretical finding is the reflectivity and the 
Kerr angle for graphite in strong magnetic fields. 



II. LANDAU LEVELS IN GRAPHITE WITH 
TRIGONAL WARPING 

Taking into account the tight-binding parameters of 
the SWMC theory, the effective Hamiltonian near the 
KH line of the Brillouin zone in graphite writes in the 
form of Rei&22. 



H(k) 



75 



vk + 71 7 4 w/c_/7o 

72 74^-/70 73v/j+/7o 

71 j 4 vk + /-fa 75 vk^ 

74«£; + /7o 73^-/70 vk + 72 



where k± = ^fik x — k y are the momentum components 
and v is the velocity parameter in the intra-layer direc- 
tions; 7j are the functions of the k z momentum in the 
main axis direction 

72 = 272 cos (2k z d ) , 75 = 275 cos {2k z d ) + A , 
7i = 2j t cos (k z d ) for i = 1, 3, 4, 

with the distance do = 3.35 A between the layers in 
graphite. The nearest-neighbor hopping integral 70 ~ 3 
eV corresponds with the in-layer inter-atomic distance 
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TABLE I: The parameters of the Hamiltonian, Eq. ((TJ), their 
values in the SWMC model, and obtained in the experimental 
works, all in meV. 
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a SWMC,Ref.-, b Mendez et al, Ref£, c Doezema et al, Ref.^ 



ao = 1.415 A and the Fermi velocity parameter v — 
1.5ao7o = 10 6 m/s. 

For the zero magnetic field, the eigenvalues of the 
Hamiltonian give four close bands. In the magnetic field 
B, the momentum projections k XiV become the operators 
obeying the commutation rule {k + , = —2ehB/c, and 
we use the relations 



k+ = y / 2\e\hB/ca, fc_ = y / 2\e\hB/ca + , 

involving the creation a + and annihilation operators a. 
We will write only one of two x, y space coordinates in- 
cluding the corresponding degeneracy proportional to the 
magnetic field in the final results. 

We search the eigenfunction of the Hamiltonian ([1]) in 
the form 



sn tp n -l(x) 



r 3 



(2) 



where <p n {%) are orthonormal Hcrmitian functions with 
the Landau numbers n > 0. One sees that every row 
in the Hamiltonian (JT|) becomes proportional to the defi- 
nite Hcrmitian function if the terms with 73 is omitting. 
We will show that the terms proportional to 73/70 can 
be considered within the perturbation theory at strong 
magnetic fields. 

Canceling the Hermitian functions from the equations, 
we obtain a system of the linear equations for the eigen- 
vector C sn 
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\ UJ^y/n — 1 



72 - e 




7i 



LU c y/n — 1 



LOi\/n 


72 - 



J 




(3) 



where the band number s — 1, 2, 3, 4 numerates the solu- 
tions at given n from the bottom, w c = vyj2\e\hB/c and 
uj 4 = 7 4 w c /7 . 

The eigenvalues of the matrix in Eq. ([3]) are easily 
found, they are shown in Fig. [T] as a function of the 
momentum k z . For each Landau number n > 2 and 
momentum k z , there are four eigenvalues e s (n) and four 
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FIG. 1: (Color online) LLs e sn for n = to 4 in four bands 
s = 1,2, 3, 4 (in dotted, solid, dashed, and dash-dotted lines, 
correspondingly) as functions of momentum k z along the KH 
line in the Brillouin zone (K = 0, H = n/2do) at the mag- 
netic field B — 7 T with the SWMC model parameters given 
in Tabl. l.The main electron transitions shown in the right 
panel below 100 meV are possible between the levels with the 
selection rule An = ±1 , see text. 



corresponding eigenvectors, Eq. ([2]), marked by the band 
subscript s. We will use the notation \sn) for levels. In 
addition, there are four levels. One of them, 



ex(n = 0) = 72 



(4) 



for n = with the eigenvector C = (0,1,0,0) as is 
evident from Eq. It intersects the Fermi level and 

belongs to the electron (hole) band near the K (H) point. 
Others three levels indicated with n = 1 and s = 1, 2, 3 
are determined by first three equations of the system ([3]) 
with Cf x = 0. The |21) level is close to the 1 10) level. In 
the region of k z , 71 / cos z 72 , where the electrons are 
located, this level has the energy 



e 2 {n = 1) = 72 



, ^c74 

'71 7o 



In the same region, the two closest bands (s 
71 > 2 are written as 



2, 3) with 



£2,3 (n) =72 - ~^(2n 
7i 7o 



7i 



1) (5) 



within accuracy of (74/7o) 2 - 

A simplest way to evaluate the corrections resulting 
from the warping 73 consists in the consideration of the 
relative Green's function having the poles at the electron 
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FIG. 2: Diagrams for the second iteration of the perturbation 
theory; corrections to the Green's function (a), corrections to 
the vertex in conductivity (b) . 



Comparing the corrections, Eq. (j9)), with the main 
contribution Eq. ([5]), we find that the perturbation the- 
ory is valid when an expansion parameter (737i/7o^c) 2 
becomes small, i.e. for the strong magnetic fields B > 
IT. The corrected 1 10} level writes 

^o) = 72+ (^) 2 £^^. (10) 

The 1 21} — level is very close to the level with n — 0, Eqs. 
@ and (HDJ). 

Our expressions for the levels with the corrections © 
and (|T0| give the same results as obtained in Ref.21 by 
the numerical method of truncating the infinite matrix. 



levels. The corrections to the levels can be found in the 
iterations 

G m+1 (x,x') = f dx"G (x,x")V(x")G m (x",x'), (6) 



where V(x) has only two matrix elements V 42 — 
w c 73a+/7o and V 24 = V 42 * in the Hamiltonian (JTJ) . The 
Green's function of the unperturbed Hamiltonian writes 



Gf(e,x,x') = Y / ^ nix) _ ) - ( 7 ) 

sn 6 £sn 

In the second iteration, we get the corrections 
d Xl dx 2 Gf{x, x 1 )V 42 {x 1 )G 22 {x 1 ,x 2 )V 24 {x 2 )Gf{x 2 ,x l ) 



and the similar term with the substitution of the super- 
scripts 2 -s-> 4. The matrix elements of the perturbation 
V are easily calculated with respect to the Hcrmitian 
functions of Eqs. (J7J, @ and we obtain for the diagram 
shown in the upper part of Fig. [2] 



To 



(e - e sn )(e - £ a ',n-a)(£ ™ e*,n) 



This correction plays an important role near the poles 
of the Green's function. Therefore, we can substitute 
e sn instead of e in the second factor of the denominator 
and represent this correction as a shift 5e sn of the poles 
(e - e sn - Sesn)' 1 with 



5e a (ri) 



^3\ 2 ^i(n-2)\C 4 n C 2 ^ 3 \ 2 



7o 



✓ I e.(»)"^(n-3) 
->2 n& i2 1 v 3 ) 



| {n + l)\C 2 sn C 4 ^V 
e s {n) -e 8 /(n + 3) J' 

where the first term should be omitted for n — 3 < 0. 
In fact, our illustration is nothing but a calculation of 
the electron self-energy and the naive expansion of the 
denominator can be indeed replaced by summarizing of 
the corresponding diagrams. 



III. CONDUCTIVITIES IN MAGNETIC FIELDS 

In the collissionless limit when the relaxation rate T is 
much less than the frequency, r«w, the conductivity is 
expressed in terms of the correlation function 

V (w) = [ dxdx'Tr {v'G ((*}+, x, x 1 ) v^G (ui-,x', x)} 



(11) 

where we should (i) summarize over Matsubara's frequen- 
cies oj m , (ii) take the trace over the Landau and band 
numbers, (iii) make an analytic continuation into real fre- 
quencies uj, and (iiii) substrate from the result its value 
at uj = (for details see Ref .22.22). 

The velocity matrices v l in Eg. (jlip are given by the 
derivative of the Hamiltonian, Eq. (JXJ) , 



dH(k) 
<9k 



(12) 



First we consider the largest velocity operators, 
Eq. (|12j) . which do not involve the parameter 73/70- 
Straightforward calculations yield two independent com- 
ponents of the dynamical conductivity 



tt/2 

ia x Juj) ir z J Ass'n (10) 

yv ' ) n,s,s' q 

x [(u + iT + A,,,,,)- 1 ±{uj + iT- A,,,,,)- 1 ] , 

where A ss > n — e sn — £y,n+i is the level spacing including 
the corrections given in Eqs. ((9]) and (fT0|). Af ss > n = 
f(e S ' n +i) — /(e S n) is the difference of the corresponding 
Fermi-Dirac functions and 

+(74/70) {c l sn c 4 n+l + c 2 sn cl n+1 ) 

is the dipole matrix element. These electron transitions 
obey the selection rule 

An = 1 , 
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and will be referenced as the strong lines. The integration 
over the Brillouin half-zone, < z < 7r/2, and the sum- 
mation over the Landau number n as well as the bands 
s, s' should be done in Eq. (|13p . The conductivity units 



Ahd 



have a simple meaning, being the graphene dynamic 
conductivity^ e 2 /4h multiplied by the number I /da of 
layers within the distance unit in the main axis-direction. 

Now we consider the terms with 73 / 70 in the velocity 
operators, Eq. (fT2"]) . Calculating the correlation function 
Eq. (fTTj) we get an additional term in the conductivity, 
which can be obtained from Eq. (fTBf with the substitu- 
tions 



1 



and with the matrix element 



d ss 'n = (73/70)^^4+2 

instead of the matrix element given by Eq. (|T4"]) . This 
transitions obey the selection rule 

An = 2 

and will be referenced as the weak lines. 

So far we considered the 73 corrections to the Green's 
function, i.e. to the levels. However, there are so-called 
vertex corrections to the self-energy shown at the bottom 
of Fig. [2j They are resulted from the quartet of the 
coupled Landau levels, which interfere while the selection 
rules An = 1 and An = 2 are allowed. For compactness, 
let us denote this quartet of given n as following 

a = \sn), b = \s , n + 1), c = \si,n + 3), d = |s l5 n + 4) . 

(15) 

where the band numbers s, s', s\, and are arbitrary. 
The corresponding corrections to conductivities write 



ir/2 



Sa xx (u,) y 
iSa xy (uj) ^ , 

* / nss sis^ q 



dz 



2^73 

7T7o 



XLj a L , b Lj c Ly d{ Ly b Lj a + °q J l^ri °c +°c L/ dJ 

x^(n + l)(n + 2)(e b - e<i) _1 (ea - £ c ) _1 
x { [(w + iT + s b - So)- 1 ±{uj + iT-e b + ej" 1 ] d ab 
+ [(cj + ir + e 6 - ec)- 1 ± (cj + «r - e b + Ec)' 1 ] d cb 

+ [(uj + ir + e d - eaY 1 ± (« + tr - e d + 6a)- 1 ] a ad 

+ [(w + ir + e d - ec)- 1 ± (w + if - e d + e c )~ l ] d cd ) 



where 



dab = [f(e a ) ~ f(e b )]/{e b - s a ) 



and f{e a ) is the Fermi-Dirac function. The terms with 
the negative radicand should be omitted while summing 
over n and all band numbers from Eq. (|15l) . 
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FIG. 3: (Color online) Real (a) and imaginary (b) parts of 
the longitudinal (xx, solid line) and Hall (xy, dashed line) dy- 
namical conductivities; Kerr angle (c) and reflectivity (d).The 
magnetic field B = 7 T, the temperature T = 0.1 meV is less 
than the level broadening T = 3.5 meV. 



IV. HALL AND LONGITUDINAL 
CONDUCTIVITIES WITH THE SWMC 
PARAMETERS 

The parameters of Eq. (jTJ) used in the calculations are 
listed in Tabl. U (see also Ref.— ). The hopping integrals 
7o t° 73 are close to the values determined in observations 
of the semiclassical ShdH effect. The Fermi energy equal 
to £f = —4.1 meV agrees at the zero magnetic held with 
the measurements of the extremal Fermi-surface cross 
sections and the masses of holes and electrons. Con- 
nections with the notation for the same parameters in 
the SWMC model are given in the "SWMC" line. The 
values of parameters 74, 75, and A determined in various 
experiments are very different, we use 75 and A obtained 
by Doezema et a£ (given in Tabl. U in the "SWMC" 
notations) and take for 74 the approaching value. In the 
quantum limit, when electrons and holes occupy only |10) 
and 1 21) levels, the Fermi energy must cross these close 
levels at the middle of the KH line. It means that the 
Fermi level becomes higher at such the magnetic fields 
taking the value ep — 1 meV. 

The results of calculations are represented in Figs. |3]|4] 
Let us emphasize that the imaginary part of the dynam- 
ical conductivity is of the order of the real part. 

One can see in Fig. [3] (a), that the longitudinal conduc- 
tivity calculated per one graphite layer tends on average 
to the graphene universal conductance. The main contri- 
bution in the sharp 16-meV line is resulted from the elec- 



5 




50 100 150 200 50 100 150 200 



Frequency (meV) 

FIG. 4: (Color online) Kerr angle and reflectivity at 10, 15, 
and 25 T. 

tron 1 21) — > 1 32) transitions (15 meV) about the K point 
where the 1 32) level coincides with the Fermi level (within 
accuracy of the width Y or temperature T). Then, the 
transitions 1 22) — > 1 21) produce the broad band. The 
low-frequency side of the band (23 meV, at the intersec- 
tion of the 1 21) level with the Fermi level) contributes 
into the 16-meV line. In the same 16-meV line, the tran- 
sitions 1 32) — > 1 33) can contribute as well if the band 1 32) 
contains the electrons. 

The next doublet at 43 meV arises from the transitions 
1 23) -> |32) and |22) -> |33) at the K point. The 68- 
meV doublet splitting results due to the electron-hole 
asymmetry from the transitions 1 24) — > 1 33) (65 meV) 
and 1 23) ->■ |34) (69 meV) at the K point of the Brillouin 
zone. 

The 89-meV line is more complicated. First, there 
are the electron transitions 1 24) — > 1 35) (89 veV) and 
|25) -> |34) (90 meV) near the K point. Besides, the 
transitions 1 1 1) — > 1 10) (95 meV) near the H point make 
a contribution as well. All these lines obeying the selec- 
tion rule An = 1 are strong. There are two weak lines 
in the frequency range. One (|24) — > 1 32) ) is seen at 55 
meV as a shoulder on the theoretical curve. Another, at 
31 meV, results from the transitions 1 10) — > 1 32) near the 
K point. 



The positions of the lines for fields in the range of 10 
- 30 T agree with observations of Refs4^. 

The optical Hall conductivity a xy (oj) in the ac regime 
is shown in Fig. [3] (a) and (b). The conductivities ct xx {uj) 
and a xy (oj) allow us to calculate the Kerr rotation and 
the reflectivity as functions of frequency [see Fig. [3] (c) 
and (d)]. It is evident that the interpretation of the Kerr 
rotation governed by the conductivity a xy (oj) is much 
more complicated in comparison with the longitudinal 
conductivity. 

The Kerr angle and reflectivity shown in Fig. 2] for 
the different magnetic fields demonstrate the strong field 
dependence of the magneto-optic phenomena. 



V. SUMMARY AND CONCLUSIONS 

In conclusions, we have evaluated the perturbation the- 
ory for the matrix Hamiltonian, which permits to cal- 
culate the corrections to the eigenvalues resulting from 
the small matrix elements particularly from the trigo- 
nal warping. We have shown that the trigonal warp- 
ing in graphite can be considered within the perturba- 
tion theory at the strong magnetic field larger than 1 
T approximately. We have found that the main opti- 
cal transitions obey the selection rule An = 1 for the 
Landau number n, however the An = 2 transitions due 
to the trigonal warping with the small probability are 
also essential. The good agreement between the calcu- 
lations and the measured Kerr rotation and reflectivity 
in graphite in the quantizing magnetic fields is achieved. 
The SWMC model parameters are used in the fit tak- 
ing into account the previous dHvA measurements and 
correcting the Fermi energy for the case of the strong 
magnetic fields. 
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